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ABSTRACT 

Using the exceptional monitoring capabilities of the MACHO project we present here 
the optical history of the LMC supersoft source (SSS) RX J0513. 9-6951, for a contin- 
uous 3 year period. Recurring low states, in which the optical brightness drops by up 
to a magnitude, are observed at quasi-regular intervals. This provides a crucial insight 
into the nature of the SSS and, in particular, a chance to investigate the poorly un- 
derstood behaviour of their accretion discs. Analysis of the high state data reveals a 
small modulation of semi-amplitude ^ 0.02 magnitudes at P=0.76278±0. 00005 days, 
a period which is consistent with the current "best" suggested spectroscopic value. 
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1 INTRODUCTION 

ROSAT observations have considerably enlarged the new 
class of high luminosity X-ray objects, the so-called "super- 
soft sources" (SSSs) characterised by their EUV tempera- 
tures (Triimper 1992). SSSs were first detected in the Large 
Magellanic Cloud in 1979-1980 with the Einstein X-ray Ob- 
servatory (Long, Helfand & Grabelsky, 1981). Until recently, 
little progress had been made in determining the exact na- 
ture of these elusive systems, the high level of X-ray ab- 
sorption rendering them undetectable in the Galactic plane. 
Most currently known SSS are therefore extragalactic and, 
as such, they are optically faint. 

The current inventory of supersoft objects (Hasinger 
1994; Cowley et al. 1996; Kahabka & Triimper 1996) is of the 
order of 11 in the Magellanic Clouds, 15 in M31 and 7 Galac- 



tic sources, with candidates existing also in MlOl, NGC253 
and M33. It appears now, through ROSAT observations, 
that the SSS do not form a strictly homogeneous class, 
some having been identified with a planetary nebula nu- 
cleus (Wang 1991), a PG 1159 star (Cowley et al. 1995) and 
symbiotic systems (e. g. Hasinger 1994). However, the bolo- 
metric luminosities of these objects (Lboi ~ 10^^ ergs~^) 
are typically an order of magnitude less than the original 
Emstem sources, CAL 83 and CAL 87. 

Among the .ROSylr discoveries are several systems, in- 
cluding RX J0513. 9-6951 (hereafter RX J0513-69) which 
show the hallmarks of the original LMC SSS. Typically, 
Lhoi ~ 10^* ergs~^ and Tbb ~ 30 eV. Spectroscopic stud- 
ies (Smale et al. 1988; PakuU et al. 1988) readily identified 
them as low mass X-ray binaries (LMXBs), yet the nature 
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of the compact object and the source of the soft emission 
proved elusive. Previously, X-rays of such low energy had 
only been observed in certain types of cataclysmic variables 
(CVs), but these accreting white dwarf binary systems are 
typically about a million times fainter than the SSS (see e. 
g. Warner 1995). 

In 1990, CAL 87 was proposed to be a black hole bi- 
nary, on the basis of a radial velocity analysis (Cowley et al. 
1990). The following year, a paper appeared advocating the 
scenario of a neutron star accrctor, shrouded in a dense co- 
coon of ionised matter (Greiner, Hasinger & Kahabka 1991). 
However, neither model seemed to provide a natural expla- 
nation for the extremely soft X-ray emission. 

The most significant progress occurred when van den 
Heuvel et al. (1992) proposed a model for the SSS which in- 
volved a white dwarf primary undergoing accretion at a rate 
> 1O"'^M0 yr-\ It was shown that steady nuclear burning 
on the white dwarf surface could produce the extremely soft 
X-ray emission at the required luminosities. However, such 
high accretion rates require a donor star more massive than 
the white dwarf, in order to sustain thermally unstable mass 
transfer. 

We have therefore undertaken a programme of optical 
spectroscopy and photometry, to test the predictions of this 
model. In particular, we focus on the LMC group, since the 
accurately known distance of these high luminosity SSSs is 
a vital key in understanding their nature. We present here 
the results of long-term optical photometry of the tran- 
sient LMC source, RX J0513-69, which was discovered in 
the ROSAT All Sky Survey (Schaeidt, Hasinger & Triimper 
1993). Remarkably, these observations were acquired as a 
serendipitous by-product of the MACHO project (Alcock et 
al. 1995a), owing to the location of RX J0513-69 in a fre- 
quently monitored field. We are thus afforded an unprece- 
dented opportunity to study the long term behaviour of this 
source which, at 1^ ~ 16—17, would normally be impossible. 



2 OBSERVATIONS 

The relative magnitude of RX J0513-69 for the period 1992 
August 22 - 1995 November 27 is presented in Fig. 1. The 
observations were made using the 1.27-m telescope at Mount 
Stromlo Observatory, Australia. A dichroic beamsplitter and 
filters provide simultaneous CCD photometry in two pass- 
bands, a 'red' band 6300 — 7600 A) and a 'blue' band 
(~ 4500 — 6300 A). However, we show here only the observa- 
tions taken through the latter filter, which is approximately 
equivalent to the Johnson V passband, since the red light 
curve is not appreciably different. The vertical dotted lines 
indicate our working definition of the 'high' and 'low' state 
sections, as used in the power spectrum analysis. 

The images were reduced with the standard MACHO 
photometry code SoDoPHOT, based on point-spread func- 
tion fitting and differential photometry relative to bright 
neighbouring stars. One-sigma error bars are shown, which 
are usually smaller than the data points (w 0.02 mag); for 
some observations with very poor seeing or transparency, 
the errors rise to ~ 0.1 mag. Further details of the instru- 
mental set-up and data processing may be found in Alcock 
et al. 1995b, Marshall et al. 1994 and Stubbs et al. 1993. 

The absolute calibration of the MACHO fields and 



transformation to standard passbands is not yet complete, 
thus the measurements are plotted differentially relative to 
the observed median. A preliminary estimate indicates that 
Am = corresponds to V « 16.3 and V — i? « 0.08, 
with systematic uncertainties around 0.2 and 0.1 mag re- 
spectively. The final calibration will greatly reduce these 
uncertainties, but they do not affect the differential mear- 
surements presented here. 

The system exhibits pronounced optical variability, 
with the most dramatic changes occurring on timescales 
of ~ 100 — 200 d. The brightness typically drops by ~ 
0.8 — 1.0 mag in only ~ 10 d, after which the light curve 
maintains a low level for ^ 30 d. It should be noted, however, 
that following the initial drop, the system usually brightens 
by 0.3—0.4 mag in the first ~ 8 d, and maintains this plateau 
level before making the rapid upward transition back to the 
high state. This pattern of variability is less obvious in the 
low state at day number ~ 1540, in which the system ap- 
pears to undergo a minor outburst, rather than exhibiting 
the step-like behaviour. 

The first and final low states of Fig. 1 exhibit a some- 
what more gradual decrease in magnitude than the others, 
taking around a week longer to drop to the faintest level. 
The overall drop is also relatively smaller by > 0.1 mag. 
The magnitude at which the system starts to enter a low 
state appears essentially constant (Am = 0) , except perhaps 
for the faint episode at day number ~ 900, which begins at 
Am = —0.2. During the high states, the brightness appears 
to show a roughly linear fading, amounting to ^ 0.2 mag in 
total, independent of the interval of time between the low 
states. We note that it is probable, both from the recur- 
rence times of the low states and from the decline of the 
light curve before day ~ 1320, that a low state was missed 
in the period ~ 1320—1370, during which there is a gap in 
the monitoring. Indeed, spectroscopic observations obtained 
in December 1992 by Rcinsch (in preparation) strongly sug- 
gest that RX J0513-69 was in an optical low state during 
this time. 



3 PERIOD ANALYSIS 

We analysed the high state sections, as defined by the verti- 
cal dotted lines in Fig. 1, for any periodic behaviour. Due to 
the fact that the intrinsic variability of RX J0513-69 of up to 
^ 0.3 mag in the 'flat' regions of the light curve can mask or- 
bital modulations, we detrended each high state section sep- 
arately, by subtracting a linear fit, before performing a power 
spectrum analysis on the combined dataset. The resulting 
Lomb-Scargle (Lomb 1976; Scargle 1982) periodogram is 
shown in Fig. 2, for a frequency space of 0.1 — 10 cycles d"'^, 
with a resolution of 0.001 cycles d"'^. We see a dominant 
peak at P = 0.76278±0.00005 d, with lesser power at 3.22 d 
and 1.45 d. However, the former of these secondary peaks 
is almost certainly a one-day alias, having a frequency of 1 
cycle d~^ less than that of the 0.76278 d period. Further- 
more, we checked the significance of the peaks by analysing 
the power spectra of randomly generated datasets, using the 
sampling intervals of the real data. Our Monte Carlo simu- 
lations reveal that the 1.45 d peak is significant at only the 
~ 1.6cr level. However, the power at P = 0.76278±0.00005 d 
corresponds to a 5cr detection, leading us to present this as 
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the true orbital period. Independent spectroscopic studies 
(Crampton et al. 1996) are consistent with this result, hav- 
ing revealed a 'best' period of P « 0.76 d. 

The high state data were folded on P = 0.76278 d to 
examine the form of the orbital modulation. Wc find that 
the data are well fitted by a sinusoid of semi-amplitude 
0.0213±0.0009 mag, suggesting a low inclination. This 
phase- averaged, folded light curve is shown in Fig. 3. We de- 
rive an ephemeris of To = JD 2448857.832(5) -^0.76278(5)£, 
where To is the time of maximum optical brightness, and E 
is an integer. 



4 DISCUSSION 

Recently, Crampton et al. (1996) reported vaxiations of only 
~ 0.3 mag, noting that the optical counterpart, identified as 
HV 5682, has historically shown variations of up to ~ 1 mag. 
Placing these observations in the context of the MACHO 
light curve, it is clear that both findings are consistent with 
the data presented here. In the more usual optical high 
states, variations of up to ~ 0.3 mag commonly occur on 
timoscales of days (see Fig. 1). Clearly, the more dramatic 
variation of ~ 1 mag can be identified with transitions to 
the rather more infrequent low states, detection of which 
requires extended monitoring of the type reported here. 

The optical luminosity in this system is expected to be 
dominated by the EUV/soft X-ray heated accretion disc, 
which has an absolute visual magnitude My ~ — 2 in the 
high state, at the high extreme of LMXBs in general (van 
Paradijs & McClintock 1995). Indeed, it is remarkable that 
all the LMC LMXBs are optically so luminous. This may 
be compared to typical values of Aly ~ +4 to +7 for CVs 
(van Paradijs 1983). Even if we apply the Warner (1987) 
empirical relation: 

Mf^^ = 5.74 - 0.259 Port(hr) (Po^ < 15 hr) (1) 

for the brightness of dwarf novae discs at maximum, using 
Porb « 18 hr, we obtain My = +1.1, substantially fainter 
than observed for RX J0513-69. This suggests that there 
is an additional source of optical luminosity in the system, 
consistent with the van den Hcuvcl et al. (1992) scenario of 
a white dwarf undergoing surface nuclear burning. 

Using our orbital period of P = 0.76278 d, we may cal- 
culate the mean density, p, of the companion star under the 
assumption that it fills its Roche lobe. Wo combine Kepler's 
third law and the Eggleton (1983) relation for a Roche-lobe 
filling star: 

a 0.6g-2/3 ln(l Q-i/3) ' 

where Rl2 is the radius of a sphere with the same volume as 
the secondary Roche lobe, a is the binary separation, and q 
is the binary mass ratio {= Mcompact/T^secondary) to obtain: 

P = p^^f^ (0-6 + <i'" + 1-'")) " g cm-, (3) 

where P is in days. For values of g < 1, as required by 
the van den Heuvel et al. 1992 model, the implied mean 
density is ~ 0.2 — 0.3 gcm~^. This is consistent with that 
of a ~ 2.5 — 3.0 Mq main sequence star (of spectral type ~ 
AO). The light from such a star would still be dominated by 



the accretion disk/compact object luminosity, since for an 
LMC distance modulus of 18.5 (e. g. Panagia et al. 1991) 
its apparent magnitude would be ~ 19, significantly fainter 
than the observed brightness of RX J0513-69 oiV ^ 17. 



5 SUMMARY 

In conclusion, we have used the remarkable monitoring ca- 
pability available as a by-product of the MACHO project to 
observe, for the first time, the long term optical behaviour 
of the recurrent SSS RX J0513-69. We see marked high and 
low states in the optical. The relation of these optical varia- 
tions to the X-ray behaviour is currently under investigation 
(Southwell et al. 1996). Indeed, ROSAT has so far discov- 
ered three X-ray on states for this transient source, making it 
the first SSS to exhibit recurrent outbursts (Hasinger 1994; 
Schaeidt, in preparation). We derive an orbital period of 
P = 0.76278 ±0.00005 d, consistent with independent spec- 
troscopic findings. Such a binary system allows solutions for 
the nature of the secondary star that are consistent with the 
van den Heuvel et al. (1992) model of unstable mass trans- 
fer onto a white dwarf. Given the small photometric orbital 
modulation (~ 0.02 mag semi-amplitude) and the extreme 
optical variability of this source, we strongly urge further 
spectroscopic observations. 
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Figure 1. The optical light curve of RX J0513-69 from 639 ob- 
servations taJjen during the MACHO project. The relative magni- 
tude is shown for the 'blue' filter, which is approximately equiv- 
alent to the Johnson V passband. Note the quasi-regular mag- 
nitude drops of ^ 1 mag. The day number is JD— 2448000. The 
vertical dotted lines indicate our working definition of the 'high' 
and 'low' state sections, as used in the power spectrum analysis. 



Figure 3. The light curve of RX J0513-69, folded on a period of 
0.76278 d. The data have been averaged into 30 phase bins, and 
are fitted with a sinusoid of amplitude 0.0213 mag. Two cycles 
are plotted for clarity. 



Figure 2. The Lomb-Scargle periodogram of RX J0513-69 MA- 
CHO time series data, excluding the low states. The strongest 
peak is at 0.76278 ± 0.00005 d. 
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